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Abstract 

This  paper  considers  the  proximity  effects  of  superconduct¬ 
ing  multilayer  films  composed  of  different  metals.  Using  the 
Silvert  and  Cooper  method,  the  relationship  between  the  critical 
temperature  of  a  superconducting  multilayer  film  with  strong 
heterogeneity  and  its  geometric  structure  is  given  under  the 
Cooper  limit. 

By  using  modern  material  preparation  techniques  such  as 
sputtering,  evaporation,  epitaxy  etc. ,  based  on  some  peoples  ' 
designs  of  atomic  scales,  we  synthesized  various  previously 
non-existent  new  types  of  materials  and  studied  their  structures 
and  properties.  This  has  been  an  important  trend  in  solid  phy¬ 
sics  in  recent  years.  For  example,  besides  semiconductor 
materials,  there  have  also  been  reports  recently  of  making  new 
materials  with  "multilayered  ultrathin  common  lattice  struct¬ 
ures"  (LUCS)  from  niobium  and  copper  possessing  different 

0 

properties  II]. 

When  compared  with  large  pieces  of  homogeneous  specimens, 
multilayered  thin  film  systems  have  two  new  factors  which  should 
be  considered.  One  is  the  appearance  of  many  hetergeneous 
interfaces.  See  Reference  [2]  for  an  analytical  discussion  of 
them.  The  second  factor  is  the  proximity  effects  caused  by  the 
many  hetergeneous  layers  being  periodically  adjacent  to  each- 
other  and  the  superconducting  electron  pairs  possessing  non¬ 
localization.  This  paper  makes  a  preliminary  investigation  of 
this  problem.  We  began  from  a  simplified  model  under 


hetergeneous  interface  conditions,  used  the  GLAG  theory  and 
gave  the  general  relationship  between  the  multilayered  film's 
critical  temperature  and  geometric  factors. 


We  conceived  that  the  multilayered  film  was  composed  from 
two  types  of  pure  superconducting  materials,  S.^  and  S2,  and 
their  thicknesses  were  separately  d1  and  d2»  In  order  to  sim¬ 
plify  the  problem  and  outstanding  focal  points,  we  assumed  that 
the  property  difference  concentrations  of  the  two  types  of  mat¬ 
erials  were  reflected  on  their  electro-acoustic  coupling 
constant  (NV) ,  separately  (NV)  and  (NV)2»  The  other  properties 
were  the  same. 

We  used  the  Silvert  and  Cooper  method  to  calculate  the  prox¬ 
imity  effects  [3] . 

The  virtual  time,  Fourier  transformation.  Green's  function 
and  abnormal  Green's  function  of  the  superconductor  satisfy  the 
following  formulas  in  the  Tc  area: 


C.(r,  r) 


-  G<?(r  -  r) 

—  ^r,driG(J>(r,  r,)A(r,)G^(r2,  r,)A*(r,)G.(r„  r), 
F«(r,  r)  -  j  Jr,Gm(r,  r,)A(r,)G'il(r,,  r'). 


CD 

(2) 


In  the  formula,  G,.  and  F  are  separately  the  completely 

/  \ 

normal  and  abnormal  functions,  is  the  zero  order  Green 

function  and  A  (r)  is  the  energy  gap  function  which  changes 
with  the  space.  They  should  satisfy  the  following  self-consist¬ 
ent  conditions: 

A(r)  -  V(r)  ■  F(r,  r),  (3) 

In  the  formula  2Lt  ^«Cr»  and  k  is  the  Boltzmann 


constant.  In  the  T  vicinity,  23 (r)  is  a  small  quantity  and 

C 

G^  and  rely  on  the  expansion  in  powers  of  23  : 

C.(r,  r  )  -  2  GL->(r.  r'),  CO 

••I 

G<T>(r,  r')  -  -  jrfr.rfrjGirCr,  r,)A(r,)G»  (r2,  r2)  •  A(rj)GL— °(r„  r'),  (5) 

r  )  -  2  Fi"^r*  r'>*  C6) 

•  *  1 

F<T>(r,  r')  -  j^.GL-^r,,  r)A(r,)G‘«  (r,  r).  (7) 

By  separately  substituting  formulas  (5)  and  (7)  into  formulas 
(4)  and  (6) ,  we  can  test  and  verify  this. 

We  can  see  from  the  above  expansion  formulas  that  fJ^  is 
the  2m-l  times  homogeneous  function  of  A  .  In  later  calcula¬ 
tions,  we  only  keep  the  first  order  terms  of  A  and  in  this  way 
formula  (3)  changes  to 

A(r)  -  V(r){tT  (  ^,GL9'(r  -  rl)A(r1)G,.#i(r  -  r,),  (8) 

m  • 

In  the  formula 

t*  ™  ( 2n  +  1  )*^i  T  *“0, 


The  zero  order  Green  function  is  known 


G'J  —  I - - - -  lexp[»^!r  -  r'  [ignw  -  —  — - —  ]. 

We  substitute  formula  (9)  into  formula  (8) 

A(r)  -  V(r)  j  drj  °'(r  -  r,)A(r,), 

In  the  formula 


(9) 


(10) 


OD 


By  assuming  each  layer  of  the  multilayered  film  is  parallel 
to  the  xy  plane,  we  can  first  carry  out  metric  integration  of 
formula  (1)  and  obtain 


:  ac*)  -  u  *)tf.  —  Z  1  a(Oe»  ( 


'2  io(z  —  z') 


In  the  formula,  Nz  is  the  fixed  region  density  of  the  z  area 
and  E.^  is  the  improper  index  integral,  defined  as 

Em(rn)  —  (  —  e*p[—  ml],  03) 

For  convenience  of  calculations,  we  make  the  following  stip¬ 
ulations:  each  layer's  coordinates  of  S1  are  kd  <  z  <  (k+a)d; 
each  layer's  coordinates  of  Sj  are  (k+a)d  <  z  <  (k+l)d,  k=0, 

+1,  +2 - ;  d=d1+d2;  a=d1/d;  JLi,  =  (VN)  2/  (VN)  . 

Each  layer  of  superconducting  multilayer  films  is  relatively 
thin.  When  we  consider  the  so-called  Cooper  limit,  that  is, 
the  thickness  of  each  layer  is  smaller  than  the  coherent  length 
of  a  corresponding  material,  at  this  time,  the  energy  gap  func¬ 
tion  of  each  level  is  approximately  a  constant,  separately  taken 
as  and  ^2*  Further,  disregarding  the  macroscopic  bound¬ 

ary  of  the  superconductor,  we  consider  the  layer  number  of  multi 
layer  film  to  be  infinite. 

By  letting  z^  and  z2  separately  be  the  two  points  in  the 
and  S2  layers  and  using  formula  (12)  for  the  two  points,  we  can 
obtain  two  linear  unrelated  equations: 


A,  -  (NV)d „A,  +  (NF^uAj. 

A,  -  rfVVWaA,  +  ftftVMj, A„ 


In  the  formula,  f^j  is  determined  by  the  precise  positions  of 
z,  and  z0.  By  taking  z,  ans  z-  in  each  period  k,  we  can  obtain 


these  two  equations.  However,  because  of  the  translational 
symmetry  of  the  system,  if  there  is  only  the  equations  between 
z i  and  z2,  then  they  are  linear  and  unrelated  and  do  not  lose 
generality.  Taking  z^  and  z2  in  the  period  of  k=0,  f^  can  be 
written  in  the  following  form: 


lm  ™  v,/2\w\.  (16) 

The  set  of  equations  (14)  and  (15)  have  non-zero  solution 
conditions  and  the  coefficient  determinate  is  zero: 


[(VX)xtn  -  U  *  l dvxyjn  -  l]  -  K™)f/»  *  fn.  O?) 

We  can  find  the  relationship  of  Tc  and  d,  a,  pc  and  (VN) 1  from 
formula  (17).  Under  common  conditions,  we  can  only  use  the 
numerical  solution  method  to  solve  formula  (17)  yet  under  strong 
heterogeneous  conditions  wherein  yU>  0,  we  can  obtain  an 
approximate  analytical  result.  At  this  time,  formula  (17) 
changes  to 

*■-(£? 

After  completing  the  integration  of  z' 

*  [«p  -  «»(**)]  +  i,  j, 

x  U  <  -  * >  -  «p(-  ^)]  ♦  S?y  [>  -  -(' 


Under  the  sign  of  integration,  the  summation  of  k  is 


/ 1. 
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(20) 


(21) 


We.  will  first  analyze  several  types  of  limit  situations 
below. 

When  a=l,  that  is,  a  homogeneous  superconductor  wherein 
the  entire  multilayer  film  is  S.^  material,  the  second  term  in 
the  curly  brackets  in  formula  (21)  is  zero  and  formula  (21) 
changes  to 

</,.)- Mb  nSn-  (:2) 

-  I  “  i 

Based  on  de  Gennes'  transformation  [4] : 

•  MbT)  £  -L-  -  In  ]^.  C3) 

I  -  I  w  I  Tr 

Formula  (18)  then  changes  into  the  T  formula  of  the  BCS  theory 


In  the  formula,  CC>D  is  the  Debye  frequency. 

When  a=0,  the  entire  multilayered  film  is  normal  metal  S 2, 
the  term  within  the  curly  brackets  of  formula  (21)  is  zero  and 
from  formulas  (14)  and  (15),  ^  A 2=0  which  is  the  expected 
result. 


When  each  layer  is  very  thin  and  d  -^0,  the  system  approaches 
being  a  homogeneous  system,  the  factor  in  the  curly  brackets 
of  formula  (21)  a  and  the  Tc  formula  changes  to 

T,  —  1.14wDexp' - ~t:  ). 

a.\  »  ' 

This  shows  that  the  system  is  equivalent  to  a  homogeneous  super¬ 
conductor  with  an  equivalent  coupling  constant  of  (aNV) . 


Under  common  conditions,  by  introducing  a  parameter  1^Q» 


formula  (21)  chancres  to 

(*K»T )  £  { 1  —  r  >  —  (Uw  > 


(26) 


ad  - 1  t3 


,  /ait  1 

/  1  —  a  '  ‘ 

1  exp(  ,  it  ) 

1  “  exp(  "  U  M 

1  —  exp 

( _  it  \ 

\  h.  ) 

(27) 


In  the  formulas,  the  selection  of  1  .  should  cause  formula 

CJ>0 

(26)  to  be  equal  to  formula  (21).  In  this  way,  the  T  formula 
can  be  written  as 


T,  *  l.Mifltjexp 


[  0VI'),(1  -  y)J' 


(25) 


By  dividing  formula  (28)  with  T  =1. 14 6>Dexp 
obtain 


we 


L 


.  T, 


1 

(NV), 


(29)  . 
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Using  copper  and  niobium  multilayered  film  as  examples.  Fig.  1 
shows  the  relationship  of  Tc/Tq  ~  a  and  x(x=d/l )  as  well  as 


When  the  curves  in  Fig.  1  are  in  x  —>00  ,  they  all  tend  to 
approach  1.  The  curve  corresponding  to  a*0.2  crosses  with  the 
longitudinal  axis  at  0.04  but  this  is  not  drawn  in  the  figure. 


Fig.  1 

In  the  above  calculations  we  did  not  consider  the  special 

effects  of  the  heterogeneous  interfaces  and  therefore  the 

results  are  basically  the  same  as  those  of  the  common  and 

simple  S-N  proximity  effects.  The  T  value  of  the  entire 

multilayer  film  system  is  between  the  large  material  T  values 

c 

of  the  two  types  of  metal.  Therefore,  if  we  find  in  the  tests 

that  the  T  value  of  the  superconducting  multilayer  film  has 

certain  differences,  for  example,  higher  or  lower  than  the  large 

material  T  of  each  constituent  element,  then  this  is  direct 
c 

evidence  of  the  influence  of  heterogeneous  interfaces  on  the 
structures  of  the  electrons. 
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